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This paper studies the way in which confinement leads to chiral symmetry breaking (CSB) through 
a gap equation. We argue that entropic effects cut off infrared singularities in the standard confining 
effective propagator 1/p 4 , which should be replaced by l/(p 2 +m 2 ) 2 for a finite mass m ~ Kf /M(0) 
[M(0) is the zero-momentum value of the running quark mass]. Extension of an old calculation of 
the author yields a specific estimate for m. This cutoff propagator shows semi-quantitatively two 
critical properties of confinement: 1) a negative contribution to the confining potential coming 
from entropic forces; 2) an infrared cutoff required by gauge invariance and CSB itself. Entropic 
effects lead to a proliferation of pion branches and a qq condensate, and contribute a negative 
term ~ —Kp/M(0) to the effective pion Hamiltonian allowing for a massless pion in the presence 
of positive kinetic energy and string energy. The resulting gap equation leads to a well-behaved 
running constituent quark mass M(p 2 ) with M 2 (0) « Kf /tt. We include one-gluon terms to get the 
correct renormalization-group ultraviolet behavior, with the improvement that the prefactor (related 
^ ' to {qq)) can be calculated from the confining solution. We discuss an integrability condition that 

guarantees the absence of IR singularities at m — in Minkowski space through use of a principal- 
is^ . part propagator. 

' PACS numbers: 11.15.Tk, ll.15.Kc 
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I. INTRODUCTION 
A. General 



i-C . The purpose of this paper is to explore possibilities of describing quark chiral symmetry breaking (CSB) with 
confining forces, in the simplest case of zero temperature and density, the only case we consider here (for finite 
temperature and density, see [l[ and references therein), without running into difficult infrared (IR) singularities. 
CSB has two essential manifestations: A running quark mass M{p 2 ) with finite M(0), and a massless pion. Any such 
description of CSB with confinement alone must also resolve the dilemma that confining forces (such as the usual 
linearly-rising potential Kpr, for string tension Kp) plus kinetic energy terms apparently have no negative terms that 
could lead to a massless pion. 

Many papers have been written on CSB with confinement, but the present approach differs from those known to 
the author. There are a number of papers 043 that make some attempt to model area-law confinement as it might 
arise in QCD (as opposed to purely phenomenological effective propagators, NJL models, and so on); all of them 
make one approximation or another, and ours is no exception. This paper differs from most of the cited papers by 
attempting to maintain covariance and avoiding the use of special gauges, such as Coulomb gauge, within the context 
of a Euclidean gap equation with confining forces. 

The present paper makes two major points. The first is to argue that entropic effects (embodied in large spacetime 
fluctuations in worldlines of pions that are composites of massless quarks) may well be a major source of negative terms 
necessary for a massless pion with area-law confinement (see Q for an early description of entropic contributions). 
These entropic effects come from the masslessness of the quarks and of the (Goldstone) pion, with the consequence 
that a pionic qq Wilson loop with a large longitudinal separation between initial and final points is highly ramified or 
branched. Anywhere along the perimeter, large transverse separation of q from q is exponentially disfavored because 
of the consequent large area-law action penalty; separation of more than about M(0) _1 practically never occurs. A 
linearly-rising potential is irrelevant for larger separations. Because the long, twisting, and narrow branches are made 
of massless constituents having a massless bound state, their action per unit length is too small to overcome entropy. 
These branches signal the formation of a (qq) condensate. [Formation of the condensate as branches of a qq Wilson 
loop may have something to do with the light-cone interpretation [9| of condensates as objects localized with respect 
to hadrons; we take no position on this possibility] On general grounds we show that the entropic effects should 
contribute a term ~ —Kp/M(0) to the pion mass, and that if only kinetic energy and linearly-rising potential terms 
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are kept, masslessness of the pion is assured if M(0) has a specific value ~ KJ . That this entropic term is negative 
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is crucial, since other negative terms, such as one-gluon exchange or hyperfine structure, may not be large enough to 
give the pion a zero mass (or in other words, to yield CSB). [We have nothing new to say about the rho meson, whose 
mass we attribute in the standard way to chromomagnetic hyperfine splitting.] 

The second point is that we can find a non-singular Euclidean gap equation for the quark running mass M(p 2 ), an 
extension of the Johnson-Baker- Willey (JBW) equation 10] to confining forces, that exemplifies both negative entropic 
contributions and an effective IR cutoff, coming from the bound on qq separation. If we use the standard effective 
propagator SirKpS^/k 4 , representing the propagator of a fictitious Abelian gluon, well-known IR singularities arise in 
a standard gap equation. This is because such a gap equation, with an open quark line, is not invariant with respect to 
Abelian gauge transformations (having nothing to do with color) of this fictitious gluon. We point out that, with the 
aid of a fictitious heavy "quark" called x> with mass M 2 ^> Kp, it is possible to construct a singularity-free Abelian 
gauge-invariant dynamics for quark CSB. This is done by using a regulated effective propagator 8ttKf S^v /(k 2 + to 2 ) 2 
in the Green's function G xq = (\T[xq{x)qx{0)}\) ■ Because G yo has a closed fermion loop it has the Abelian gauge 
invariance that leads to a cancellation of the IR singularities @, [ljj, so the to = limit exists, at least as long as there 
is a mass gap, that is, that there is CSB for the quark. In such a case gauge invariance gives a natural IR cutoff in 
the dynamics at a momentum scale ~ M(0). 

Part of this cancellation of singularities was given long ago [3] in another closed-loop process, the pion Bethe- 
Salpeter equation. There it was shown that the limit m —> in the sum 2M + V(r) of the on-shell running quark 
mass and the static potential based on the regulated effective propagator exists and is finite; the potential is just Kpr. 
However, any finite terms surviving at m = were not investigated. In the present work we show that there is such 
a finite and negative term in the mass 2M of order —Kp/M(0) that we now identify with an entropic contribution. 
Such a negative term is essential for producing a zero-mass pion. In the mass-regulated potential considered by itself 
there is, at r — 0, a negative term —Kf/ttl. We see that in the m = limit, the regulator mass reappears as a term 
that we can interpret as entropic for a finite value of m ~ M(0). Moreover, finite m leads to an IR cutoff of the 
type we expect, which means that there is no reason to insist on a linearly- rising potential much beyond distances 
~ M(0) -1 . The result is a non-singular description of CSB using the mass-regulated potential with finite to, that 
yields M 2 (0) » Kf/tt. With only confining forces, the UV asymptotic behavior is M(p 2 ) ~ 1/p 4 , but we add one- 
gluon effects to get the known [12J renormalization-group (RG) behavior, with the additional feature that we can 
calculate the prefactor from the confining effects. 

As a subsidiary point, we investigate the contributions to CSB from a one-dressed-gluon JBW equation. If, as 
previously argued from this equation [131 ] . gluonic effects are not enough for CSB, they are still important because 
they enhance the zero-momentum quark mass M(0) and give the dominant effects in the ultraviolet (UV). 

Finally, we briefly mention another approach to removing the singularities of confining forces, by using a principal- 
part propagator in Minkowski space. 



B. The effective confining propagator and the gap equation 

The effective Euclidean propagator (having nothing to do with the true QCD gluon propagator!) that we will use 

is: 

D eff {k),^8, v D eff {k)- D eff (k)= ( ^ a) , . (1) 

with a finite value of to. The to = limit, or 8irKF/k 4 , of this effective propagator mimics an area law when used in 
a world-line action: 



dx,j, j>dy v SnvDeff(x - y;m = 0) (2) 

where the integrals are over closed world lines and D indicates a Fourier transform. It is precisely an area law in 
d = 2, where the (massless gluon) propagator is ~ S^u/k 2 , and it yields approximately an area law in d — 4. (For a 
discussion of how this approximates an area law in d = 4 see [Hj].) It is gauge-invariant in the sense that if derivative 
terms are added to the propagator they contribute nothing to the action, because of the integration over closed loops. 
For finite m this gauge-invariant action does, of course, depend on to, but for infinitesimal m a small change is an 
Abelian gauge transformation Q and this closed-loop action is both independent of to and non-singular at m = 0. 
This choice of massive effective propagator is not unique; as an example, one might try the two-mass propagator: 



D eff (k),^S^D eff (k); D eff (k)= {k , +ml){ ^ +ml y (3) 
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If both mi, m% are close to M(0) there is little difference from Eq. (flj. But if either is very small or very large compared 
to M(0), either unphysical features appear or CSB is absent. For example, at mi = the two-mass propagator is 
effectively an attractive Coulomb potential with a very strong coupling that has no physical basis. So we use only the 
propagator of Eq. ([Tj in what follows. 

The Euclidean JBW equation we use to describe confining effects is: 

where D e f f(p — k) is the effective propagator of Eq. ([TJ . 

It is not hard to see that, just as for the gluonic gap equation 15], the confining gap equation for M(p 2 ) is closely 
related to the gauge-invariant Bethe-Salpeter equation for the pion, with confining forces. Recall that the usual 
derivation of the JBW equation, for gluonic exchange rather than confinement, comes from the Schwinger-Dyson 
equation, with zero bare (current) mass, for the CSB mass term M(p 2 ) in the quark propagator. This equation is 
coupled to the equation for the coefficient of f in this propagator and to various vertices related to the propagator by 
Ward identities, and in its full glory the coupled equations are quite complex (see, for example, a very recent study 
of the full Schwinger-Dyson equation for the Landau-gauge quark propagator 16]). Extensive use of these Ward 
identities and other tools show [1 51 ] that satisfaction of the full quark self-energy Schwinger-Dyson equation implies 
the existence of a zero- mass pion, whose wave function is closely related to the mass function M(p 2 ). The confining 
JBW equation we use is similar, but with gluonic lines replaced by the effective propagator of Eq. ([TJ). For simplicity, 
and because further accuracy is not warranted, we consider here (as many authors do) only a simplified form of the 
self-energy equation in which the coefficient of p is unity. As with any effective propagator, we only use it at one-loop 
level. 



C. Gluon exchange in the JBW equation 

In addition to these studies of confinement we briefly explore the usual Euclidean JBW equation for one-gluon 
exchange, but with a massive gluon [131 ] . One reason to do this is to set the stage for techniques used in the confining 
JBW equation. Another is that within the general framework of the gap equation we use for confinement, but with a 
one-gluon propagator and running charge rather than the effective confining propagator, it appears that one-dressed- 
gluon exchange does not yield CSB for quarks [Hj], but would do so for adjoint fermions (coupled more strongly by 
a factor of 9/4 in QCD), as shown in lattice simulations [l7| . This one-gluon result is based on the non-perturbative 
generation of a dynamical gluon mass [HI, [l!| m g of about 2A where A « 300 MeV is the QCD mass scale. The 
main effect of the dynamical mass is to reduce the zero-momentum strong coupling a s (0) to about 0.4-0.5 (with no 
quarks). The conclusion that one-gluon CSB does not give CSB is far from unasssailable, because of approximations 
in the gap equation itself and possible inaccuracies in the claimed dynamical mass and running charge. Indeed, a 
preprint that came out as this paper was being written up [161 ] , based on an extensive and complicated study of the 
gap equation in Landau gauge and using Landau- gauge form factors from lattice simulations as input, claims that 
inclusion of rather subtle ghost effects does lead to gluonic CSB without confinement. As these authors acknowledge, 
the CSB mechanism they find is not very strong, and they also confirm that CSB is absent within the approximations 
of Ref. [l3T ]. which has no explicit quark-gluon vertex corrections or ghost contributions. 

In our opinion the fate of gluonic CSB remains to be determined. There are lattice simulations [201 ] claiming that 
in 5/7(2) lattice gauge theory, confinement by center vortices is both necessary and sufficient for CSB, suggesting 
that standard gluonic effects are not the CSB driving mechanism. There are also simulations (for example, |21| ) 
showing that the CSB phase transition temperature is quite close (not necessarily identical) to the deconfinement 
phase transition temperature. Another recent preprint [22j muddies the waters further, claiming that in SU(3) lattice- 
gauge simulations removal of center vortices removes confinement but not CSB, contradicting the SU(2) simulation 
results of |2(|. We do not know how this puzzle will be resolved. In order to be consistent we continue to use the 
one-gluon gap equation of Ref. [l3T ] , and add it to the confining gap equation to recover known [l2l [23| renormalization 
group (RG) and operator-product results in the deep UV, which are that M(p 2 ) ~ const. [In p 2 ] a ~ 1 /p 2 (here a is the 
Lane [12j constant) for a constant undetermined by the RG or OPE that is proportional to (qq) . In our combination 
of confinement and one-gluon effects we can estimate the unknown constant. 

D. Another path: The integrability condition in Minkowski space 

There is another way of treating to as a regulator, at least in Minkowski space. This is to define the effective 
propagator analogous to Eq. Q as a principal-part propagator, which has two important effects. The first is that it 



FIG. 1: Two-loop graphs whose sum has Abelian gauge invariance for the confining propagator. The thick line is the x field; 
the thin line is a quark with zero current mass; the curly line is the "area-law" propagator. 



makes the quark action of Eq. (|2j real for any m; there is no good physical interpretation of an imaginary part of this 
action. The second effect is that use of a principal-part effective propagator leads to an integrability condition that 
allows for the limit m — 0. We save a detailed investigation of this propagator, and its extension to Euclidean space, 
for a later paper. However, it is important to point out that the principal-part effective propagator obeys a simple 
condition, called the integrability condition, which is that the spacetime integral of the confining effective propagator 



is a finite dimensionless constant independent of any regulation scheme used to define the effective propagator. This 
real propagator is half the sum of the advanced and the retarded propagators, or the Wheeler-Feynman propagator. 

Our major claim is that whether one uses technique allowing m to be thought of as a true regulator, or whether 
one identifies a —Kp/m term as a finite entropic contribution in the sort of gap equation we consider here, proper 
choice of m will lead to correct CSB predictions from such a gap equation. 



Section HTl shows how to extract the JBW equation from the Green's function G xq . Section Hill develops the entropic 
argument, basically that configurations of a Wilson loop with massless quarks heavily favor highly-ramified structures 
with pion branches, rather than configurations where the q and q can separate by large distances (which would be 
allowed, in principle, in the quenched approximation where confining strings cannot break). Section llVI calculates 
a previously-neglected regulator-independent term in an old argument @ concerning the cancellation between a 
regulator in a 1/fc 4 potential and this regulator as it appears in the on-shell quark mass M = M(p 2 — M 2 ) calculated 
from the one-loop Feynman graph with this potential. This term, which is negative and ~ —Kp/M, is crucial to 
find a massless pion given that naive confinement contains only the positive area-law potential and kinetic energy; 
we identify it with entropic contributions. (Negative terms such as one-gluon exchange and hyperfine splitting are 
too small, or else they would give CSB without confinement.) Sections El IVII and [ VIII discuss Euclidean JBW gap 
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equations for, respectively, massive one-gluon exchange; a confining potential with finite regulator mass m ~ K F ; 
and a JBW equation for both confining forces and for one-gluon exchange. Section [Villi introduces the integrability 
condition for removing singularities of a confining effective propagator in Minkowski space. 



We argue that it is possible to extract a JBW gap equation from an Abelian gauge-invariant amplitude that is non- 
singular and that has all the physical effects (entropy, IR cutoff) represented with reasonable quantitative accuracy. 
This extraction is in somewhat the spirit of the pinch technique [la, Il9j. in which (color-) gauge-invariant Green's 
functions are extracted from some gauge-invariant object such as the S-matrix. Consider the gauge-invariant heavy- 
light Green's function G xq = (\T(xq(0)qx( x ))\) > as shown in Fig. [TJ As appropriate for an effective propagator, we 
keep only graphs with a single effective-propagator loop. One can show either by direct calculation of graphs or by 
using the techniques of jTIJ that such a gauge-invariant color-singlet amplitude is of the form: 




(5) 



E. Organization of the paper 



II. ABELIAN GAUGE INVARIANCE, INFRARED CUTOFFS, AND ENTROPY 




(6) 
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Whatever the value of m, T^ u is gauge invariant and conserved: k^T^ = 0. A mass gap plus conservation imply that 
Tp V must vanish at k = 0, and so the m — singularities cancel in the sum of all potentials and masses. The JBW 
gap equation comes from dropping all terms referring to the held \ m the sum °f graphs, after cancellation of gauge 
artifacts. 

Note that in practice if one uses the gauge-invariant description of CSB through the Green's function G xq little 
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is lost by keeping a finite m, provided that m < K F , simply because smaller momenta are cut off. This is exactly 
what we do in finding the JBW gap equation, because a hnite m ~ M(0) incorporates both the gauge-invariance IR 
cutoff and the entropic term. 

We see these properties in the static potential of this effective propagator, which is: 

V(r) = = + - e— ) (7) 

m mm 

where the first term on the right of the last term is V(r = 0). We will identify this negative term with entropic 
effects, choosing a specific finite value of m. Note that the small-mr limit of the remaining finite term is the usual 
linearly-rising potential Kpr. The r dependence at distances >~ 1/M(0), which in this case is virtually flat, should 
not matter much because separation of the q and the q beyond this distance is very improbable. 



III. ENTROPY AND A qq CONDENSATE 



It is rather easy to understand the basic properties of area-law dynamics for heavy quarks \i where heavy means 
that the quark mass M x obeys M 2 ^> Kp. (We assume these quarks are quenched and that there are no other 
matter fields.) Such quarks move in essentially classical paths, nearly straight lines, even for times T that can be 
large compared to the QCD time scale, so that the spatial separation R of a XX P a i r can be specified in advance 
with small changes coming from the quark dynamics. An area law simply means that the expectation value of the 
RxT Wilson loop describing this configuration is about exp[— KpRT]. There is no question of a condensate of these 
quarks, since their paths (in Minkowski space) have essentially no backward-moving segments, necessary for a (xx) 
condensate. Moreover, the entropy of the flux sheets confining these quarks is not large compared to the action (or 
else the 't Hooft criterion would be violated, yielding both confinement and dual confinement). 

For light quarks q, with (current) mass obeying M 2 <C Kp, things can be very different. (Again assume the 
quenched case, with no other matter fields; string breaking is impossible.) Suppose we fix the time T between the 
initial and final qq configurations, each for simplicity taken to be at the spatial origin. What kind of paths occur in 
the path integral for a quantity such as (T[q'j5q(0)q'y^q(T, 0)]) , and what kinds of areas do these paths have? 

The answer is that the paths are highly ramified (branched), as shown in Fig. [3] (For clarity we do not show 
the fine-scale movement of the quarks about each other, which is sketched in Fig. [3]) We can no longer specify the 
average spatial separation R of these light quarks, which must be calculated. An average Wilson loop resembles 
a highly-branched polymer, with most of the branches representing pions with mass proportional to M q . Note that 
going backwards in time is not at all hindered, so that there is no bar to forming a non-zero value of (qq). A typical 
configuration will have an overall length L ~ T 2 /£ , but any particular branch will have a separation between the q 
line and the q line of 0(£). Here the correlation length I is related to the CSB-generated constituent quark mass M 
by Kp £ ~ M (see, for example, an elementary discussion in [24|). So the area of the Wilson loop is 0(L£). In terms 
of the overall length L, this looks like a perimeter law, although the actual area can be large because of ramification. 

The upshot of this discussion is that when Kp ^> M q , pionic Wilson- loop configurations that look sheet- like (that 

is, for some qq separation R that is large compared to K F 1 ) are highly suppressed. This is because the available 
configurational entropy of the flux sheets cannot overcome the action penalty in the exponentially-small area law 
contribution to the Wilson loop VEV. But because the action penalty for forming pions from light quarks is small, 
entropy can dominate. Roughly speaking, the entropic term to be added to the area-law term gives a result of the 
form: 

(W) ~exp[-i. F T 2 + ^ln(2d-l)] (8) 

in spacetime dimension d. Since the overall length L scales as T 2 /£ this can also be written: 

(W) - exp[-K F L£+ - hx(2d- 1)]. (9) 

(We use a standard approximation for the entropy, as counting the ways links of length £ can extend themselves on a 
hypercubic lattice with no backtracking. There are other terms in the action coming from quark-pion vertex effects, 
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(0,0) 

FIG. 2: A schematic of a Wilson loop for a pion characteristic of light-quark dynamics; the thick line symbolizes the original 
qq loop, and the thin lines symbolize pions. Because the loops are narrow we do not show their individual q and q lines. A 75 
is understood at each end point and every three-vertex. 



and other terms in the entropy coming from counting the ways the Wilson loop is ramified; we do not discuss them 
here.) With I ~ M/Kp the entropic term contributes a term ~ —Kp/M to the action density, or energy. The Wilson 
loop will ramify until the two terms are approaching equality, at which point other physical effects take over. 

Even though rho mesons (for example) are made of virtually massless quarks, they are heavy in part because of 
QCD hypcrfine interactions to the point that their fluctuation entropy cannot overcome their action. A standard 
estimate of hyperfine splitting is: 



M p -M^^P-\m\ 2 (10) 
which is about 700 MeV for parameters that we use in this paper. 



IV. HOW CAN A PION BE MASSLESS WITH CONFINING FORCES? 



For purposes of the following heuristic discussion, we assume that one-gluon effects are too weak to produce CSB 
by themselves, and omit writing them. 

A long-standing problem of CSB via confinement is that a linearly-rising potential such as Kpr is positive; when 
added to the positive kinetic energy of the quarks, how can one get a zero pion mass bound state? The answer is the 
negative entropic contributions of the last section. We saw there that the effective confining propagator of Eq. ([1]) 
has such a term, although there was no explicit reference to entropy. 

There is another way to find such a negative term. The Minkowski-space version of the gap equation, Eq. 
and its IR divergences were studied in Q, but without mentioning an important point that we take up here. In the 
right-hand side of this Minkowski-space gap equation set M(k 2 ) to a constant value M, and calculate the integral 
on-shell, that is, at p 2 = M 2 . The result is the running mass on-shell, or M(p 2 = M 2 ). In a gap equation this should 
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FIG. 3: A schematic closeup of the qq Wilson loop in Minkowski space, showing the propagation of massless particles along 
their light cones under the influence of a confining force (after Very massive quarks would only propagate in the forward 
timelike direction, but light quarks can have spacelike legs, signifying formation of a condensate. 



be the same as the input, which leads to, as given in 0]: 



M 



2m 



(11) 



where the omitted terms either vanish at m = or arc independent of it. The static potential is just that of Eq. (|7|), 
and for small to it looks like: 



V(r) 



-— + K F r + 0(m). 
m 



(12) 



So 2M + V(r) is free of infrared divergences in the to = limit. What was not explored in this earlier paper was the 
to- independent term of Eq. (|lip . We have now calculated this term; adding it to 2M one finds: 



2M + V(r) = K F r - 



3K F 



(13) 
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In effect, although the —Kp/m term in V(r) cancels, it reappears in finite form as —3Kf/(ttM). At first it may seem 
odd that the mass operator gives a negative term to the sum 2M + V(r), but something must do so, or it will be 
impossible to find a zero- mass pion (in the present crude approximation). 
Look at the relativistic pseudo-Hamiltonian 

H= P + K F r-^. (14) 
Substitute p — > 1/r and minimize on r to find a variational approximation 
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There is a zero- mass bound state when M = 3K F /(2tt), or about 220 McV. Of course the estimate coming from 
Eq. (|15[) is only qualitative, and in the real world there are other negative terms to be included, including gluon 
exchange and hyperfine structure, but these are not the dominant negative contributions. (If they were, the one-gluon 
JBW equation would have yielded CSB, but it is perhaps plausible that it does not.) 
If we identify the two negative constants in Eqs. (|13I12I) . we find that m is close to M: 

nM eta 
m=— . (16) 

This is not necessarily an accurate value, so in the Euclidean confining gap equation using the Euclidean effective 
propagator of Eq. ([I} we set m — aM for a range of a s» 1. But first, both to set the stage for a Euclidean 
phenomenology of confinement CSB and to illustrate the argument that one-gluon exchange may not be enough for 
quark CSB, we briefly review the arguments of [l3j concerning one-gluon exchange. 



V. MASSIVE ONE-GLUON EXCHANGE 



There is no doubt that the JBW equation for QCD with a massless gluon is correct and useful for asymptotically- 
large momenta [T^], and in fact we will recover these results as an addition to CSB by a p urely confining force. 
However, it can be questioned whether one-gluon effects in the IR can produce CSB for quarks |13j . 

Long ago the author argued that the infrared singularities of QCD, coming from asysmptotic freedom, had to be 
cured by the generation of a dynamical gluon mass [18| . The pinch technique (PT) and the gauge technique were 
used to enforce gauge invariance in off-shell Green's functions for non-Abelian gauge theories. (For a comprehensive 
treatment of the pinch technique see Qjl; the gauge technique constructs vertices as functional of self-energies such 
that the Ward identities are exactly satisfied. These vertices are approximate, but expected to be asymptotically 
accurate at small momenta.) A one-dressed- loop PT approximation showed that "wrong-sign" asymptotic freedom 
problems were cured by generation of a dynamical gluon mass of about 600-700 MeV or so. In recent years such a 
dynamical mass has been abundantly confirmed by lattice simulations and more sophisticated PT treatments; see 
[l9l | . As for the running charge, [18J gives the following approximation for Euclidean (spacelike) momenta: 

S ^ ^ = 61n[(fc 2 +4m2)/A 2 ]' "'^ = 47r&ln(4m 2 ,/A 2 ) ^ 

where m g is the dynamical gluon mass, A the QCD mass scale, and, for gauge group SU{N) with Nf flavors, b is the 
one-loop coefficient in the beta-function: 

UN-2N f , , 

b =^8^- (18) 

(We use this running charge only only for spacelike momenta, so the singularity for timelike momentum is irrelevant. 
There is a more complicated modified form [25j | that is free of singularities in the timelike regime, and it agrees 
rather well with the above form for spacelike momenta.) In this paper we use, consistent with lattice determinations, 
phenomenology, and more sophisticated treatments of the PT Schwinger-Dyson equations 119], m g = 2 A s» 600 MeV. 

Let us accept these PT results, although the final word has yet to be said on their quantitative accuracy, and ask 
what they have to say about one-gluon CSB. Recently the author argued [13[ that the gluon mass and relatively small 
coupling were, in fact, too small for quark CSB via the gluon-exchange JBW mechanism, and confinement had to be 
the main source of CSB, as we now review. The linearized massive-gluon JBW equation is, in Landau gauge: 

M(p)= (2^y dJb [(p-ib)=»+ m g]^ (19) 
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where C2 is the quark Casimir eigenvalue (4/3 in 5(7(3)) and m g the gluon mass. It turns out that accounting for 
a gluon mass does two things: 1) It makes the gap equation finite at zero momentum; 2) and more important, it 
bounds the IR running coupling. If the mass is too small, the running charge gets unacceptably large, as judged by 
phenomenology and solutions to the PT Schwinger-Dyson equations. 

In the IR we can replace the running charge g 2 (q 2 ) by its zero- momentum value, which we call simply g 2 , and then 
integrate over angles, with the result: 



dfi t 



(p — k) 2 + m 



= K(k;p) = 



4tt 2 



p 2 + k 2 + m 2 + [(p 2 + k 2 + m 2 ) 2 - Ap 2 k 2 } 1 / 2 ' 



(20) 



However, this kernel does not yield a simple differential equation. Earlier, it was proposed [1< 
the angular integral by: 



K(k;p) « 2vr 2 



9(p 2 -k 2 ) 6(k 2 -p 2 ) 



k 2 



mi 



K(k,p). 



|25( to approximate 



(21) 



Numerically the approximate kernel K is, on the average, about 20-30% larger than the true kernel for IR momenta 
(< m g ), but it approaches the true kernel in the UV. We ignore this IR discrepancy, because it is in the direction to 
reinforce our conclusion that one-gluon exchange is too weak for CSB, and because the primary use of the one-gluon 
JBW equation will be for large momenta. Using the approximate kernel K and the appropriate arguments for the 
running charge yields an integral equation: 



M(p 2 ) 



iC 2 g 2 



dk 2 



k 2 M(k 2 ) 
M 2 (k 2 



16tt 2 ./ k 2 
There is a corresponding differential equation: 

2M'(p 2 ) \M{p 2 ) 



6{p 2 



(k «-> p) 



(22) 



M'V) 



(V 2 



3C 2 g 2 
16tt 2 



This is nothing but the original JBW equation, with the variable p 2 



M{p 2 ) = const, (p 2 + m 2 ) 



2\v ± 



= 0; A = -r^r. (23) 
in place of p 2 . It has power- law solutions 

(24) 



1 



v ± = ~{-l±[l 



4A] 



l/2i 



If the zero-momentum coupling is too small, there is no CSB. The standard analysis is that the critical coupling is 
the point at which the square root in Eq. (|24p turns imaginary, and that there is CSB for couplings larger than this 
critical value. Then one-gluon CSB, in this approximation, requires: 



a s (0) 



g 7r 
4^ - 3CV 



(25) 



With C2 = 4/3 this yields a s (0) > 0.8, approximately, somewhat greater than the value 0.5 given by Eq. (jT7|) . Taking 
account of the difference between the true massive kernel and the approximate kernel would change the critical value 
of a)s(0) to about one. It seems likely, then, that one-gluon exchange is not strong enough to drive CSB. 

There are two results of this subsection: First, the suggestion that ordinary one-gluon exchange is too small to drive 
CSB for quarks (but likely large enough for adjoing fermions, where the critical coupling is only 4/9 as large, because 
of the Casimir eigenvalue |13l|). Second, one can derive a regulated gap equation for an area law by differentiating 
with respect to m 2 , and replacing the coupling g 2 by —8itKf- We will pursue this second course in order to propose 
another form of an area law potential, partly based on the considerations of Sec. IVIIIl 



VI. THE CONFINING GAP EQUATION 



The issue is to solve the equation: 
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Clearly, we can find the S-wave projection of the gap equation by differentiating the massive kernel of Eq. ([20)) 
with respect to the gluon mass. But (as with the one-gluon equation) the resulting gap equation can only be studied 
numerically. Again without losing any essential accuracy, we will instead differentiate the approximate massive one- 
gluon S-wave kernel K of Eq. (|21[) with respect to m g , replace the gluon mass m g by m, and make some other obvious 
changes. This yields a Euclidean gap equation: 

2K F f p \ l2 k 2 M(k 2 ) 



Tr{p 2 + m 2 ) 2 J dk \ 2 + M 2 (k 2 ) 

2 \ I J (J2\ 



2K F 



dk 2 



k 2 M{k 2 ) 



n 2 + m 2 ) 2 
= J>(P 2 ) + J<(P Z ) 

where J>(p 2 ) is the integral from to p 2 . The corresponding differential equation is 



(fc 2 +m 2 ) 2 (k 2 + M 2 {k 2 )) 



(27) 



2\ll 



3M(p 2 )' AK F 



p 2 M(p 2 



= 0. 



(28) 



p A + to- 2 7r \_{p 2 + m 2 ) 3 [p 2 + M 2 (p 2 )] 
This is really a family of differential equations, as we see by writing them in non-dimensional form. Write 

_ 2n,r2\ ;i/r — ntfjz _ Q) ; (29) 



M(p 2 ) 
m 



Mf{u = p 2 /M 2 ), M = M{p 2 
aM 



to find: 



/"(«) + 



3/'(u) , 4K F 



u + a 2 



ttM 2 



*/(«) 



a 



2\3l 



0. 



(30) 



The boundary conditions are /(0) = 1, /'(0) = 0. From Eq. (ITB1) we expect awl. Suppose for the sake of argument 
that a is fixed even if M changes; this is reasonable, given that m is not an externally-imposed quantity but one 
which actually should scale with M. Then note that the dimensionless coupling parameter K F /M 2 appearing in 
the differential equation is not known if only the differential equation is available. It has to be determined from the 
integral equation at zero momentum: 



M 



2K F 



du- 



uf(u) 



2K F 



(u + a 2 



f 2 (u)] 



1(a). 



(31) 



Either the integral or the differential equation yields a large-momentum falloff M(p 2 ) ~ 1/p 4 . That the UV falloff 
is faster than what is expected from the OPE and the renormalization group is not a problem; the required behavior 
M(p 2 ) ~ l((lnp 2 ) a /p 2 [12] follows from one-gluon exchange, which we take up in the next section. While there is 
no need for an area law to give the actual UV behavior, it is necessary that the area-law UV behavior be no slower 
than prescribed by the OPE and the renormalization group. The correct RG behavior will be reinstated by adding 
the one-gluon terms of Eq. (|T9")) . 

For small momentum one finds: 



/(«) = [1 



2K F 
3ttM' 



-hr) + 



(32) 



showing that the running mass changes but little near zero momentum. 

It remains to determine A/, a quantity that comes entirely from J<. Given the family of solutions to the differential 
equation, we can now estimate the mass by imposing the integral condition of Eq. pip . It turns out numerically that 
M does not change very rapidly with m in the vicinity of m — M . Calculations for the range 0.8 < a < 1 yield mass 
values in the range M 2 — (0.6 — l)K F /ir, with smaller a corresponding to larger M. (The limit a or m = gives 
M diverging like ln(l/m)). We show the case a = 0.9, for which M w 0.9y/K F /n m 230 MeV, in Fig. 0] Note that 
p 2 /M 2 fa 15 corresponds to p 2 s» 1 GeV 2 . 

As we have seen, at fixed a the only dimensionless strength parameter in this zero-temperature, zero-density problem 
is K F /M 2 . Unlike the dimensionless coupling constant of one-gluon exchange, this parameter cannot be tuned, since 
M ultimately is determined by an integral equation. So there can be no phase transition between CSB and not-CSB 
at some finite value of this parameter. If one tries to find a phase with no CSB, so that M somehow approaches zero, 
the dimensionless strength parameter becomes infinitely strong, and there surely must be CSB. 
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5 10 15 20 

FIG. 4: Numerical calculation of the running quark mass f(u) vs. u, for M 2 = 0.8K f /n and a = 0.9. 



VII. A GAP EQUATION WITH BOTH AREA LAW AND ONE-GLUON TERMS 



One-gluon exchange determines the UV asymptotic behavior through the running charge, so we reinstate g 2 instead 
of the zero-momentum value g 2 in the one-gluon equation. Then the area-law plus one-gluon integral equation is: 



M(p 2 ) 



dk 2 



dk' 



k 2 M(k 2 



2K P 



3C 2 g(p 2 ) 



2Y2 



k 2 + M 2 (k 2 ) 

k 2 M{k 2 ) 
{k 2 + M 2 (k 2 )) 



ir(p 2 



+ M 2 ) 2 
2K F 



16tt 2 (p 2 + m 2 ) 
3C 2 g(k 2 ) 2 



(33) 



7T(fc 2 + M 2 ) 2 



167r 2 (fc 2 + m 2 ) 



J>(P 2 ) + J<(P 2 ) + Kyip 2 ) + K^p 2 ) 



where J>, J< are defined as before (see Eq. (|2"7jl ) and K > ,K < refer respectively to the integrals from to p 2 and 
from p 2 to infinity of the one-gluon kernel. Note that for QCD the one-gluon term begins to dominate the confining 
term at a momentum p 2 of order 2Kp /a s (0). 

One-gluon corrections to confinement are of two types. The first correction is in the IR; at zero momentum it 
comes entirely from K < . Let us now call the mass coming solely from confining effects M c . The correction to M c is 
approximately: 



M = M r 



1 



[1 - all] : 



h = 



du 



uf(u) 



(ii + 7)(w + / 2 (M))ln[/3( M + 4 7 )] 



where j3 = M 2 /A 2 , 7 = m 2 /M 2 , and a is the Lane constant: 



9C 2 



a = 



3C, 



16tt 2 6 UN-2Nf' 



(34) 



(35) 



As expected, one-gluon corrections increase M, since one-gluon effects work in the direction of producing CSB. For 
the case shown in Fig- HI plus f3 « 0.7, (3j = 4, we find 1\ w 0.22, and with (for no quarks, SU(3)) a = 4/11 the mass 
is increased by a factor of 1.1, approximately, or to about ^Kf/tt ~ 250 MeV. 

The second, and more important to us, is the UV correction, where the one-gluon term dominates. There is one 
useful simplification: The term K < (last term on the right of Eq. ([33])) is non- leading by one power of \np 2 in the 
UV, and we will drop it. Define M c (p 2 ) as the solution to the pure confinement equation (12"?| with a kernel J c (p; k): 



M c {p 2 



J c {p;k)M c (k 2 ); J c {p;k) 



2K F k 2 



6(p 2 



{p 2 +m 2 ) 2 [k 2 + M 2 (k 2 )Y 



(k O p) 



where 



dk 2 



(36) 



(37) 
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Write the solution to Eq. (|33)l . without the K< term, as M(p 2 ) = M c (p 2 ) + Q(p 2 ). We are only interested in the UV 
behavior of this equation, which then can be linearized in Q, and takes the form: 

M c (p 2 ) + Q(p 2 ) = Jj c (p; k)[M c (k 2 ) + Q(k 2 )} + J K> (p; k)[M c (k 2 ) + Q(k 2 )} (38) 

where the running mass in the denominator of K > is M c (k ): 



K>{p;k) 



ak 2 9(p 2 - k 2 ) 



p 2 



H^^W + M 2 c (k 2 )] 



(39) 



A 2 

In Eq. ([38]) the M c on the left cancels against the J C M C term on the right, leaving: 

Q = J C Q + K > (M C + Q) (40) 

using a streamlined matrix notation, with momentum arguments and integral sign suppressed (this should not lead 
to confusion of M c as we now use it with the mass as defined in Eq. ([3710 . We will solve this in the UV as a power 
series in K > . It will turn out that Q(p 2 ) ~ (lnp 2 ) a ~ 1 )/p 2 in the UV, and it is straightforward to see that the term 
J C Q vanishes more rapidly than this, by a power of p 2 . So we drop the J C Q term in Eq. (|40|) . This leaves: 

Q = K>M C + K>Q = - l — K>M C = (1 + K> + (iV>) 2 + . . .)K>M C . (41) 
1 - K > 

Defining the inverse (1 — -fT>) _1 is slightly subtle, because the integral of K > over a function depends on how 
rapidly the function vanishes in the UV. In particular, the confining solution M c vanishes like l/p A in the UV, which 
means that the function K > M C vanishes like l/p 2 \np 2 : 

K>M c (p 2 ) = ° [ P dk 2 k 2 M c (k 2 )[k 2 + M 2 (k 2 )}-\ (42) 

[p 2 + m 2 g ] In [^-£3-*] Jo 

Because of the rapid vanishing of Mc{k 2 ) in the UV we can the upper limit in the integral from p 2 to infinity. Then: 

K> M c (p 2 ) ee a/c rp2+4m2l (43) 

with 

poo 

I c = dk 2 k 2 M c (k 2 )[k 2 + M 2 (k 2 )}- 1 . (44) 
Ja 

But for any function F(p 2 ) that behaves like 1/p 2 hip 2 at infinity the integral K > F behaves like (lnlnp 2 )/p 2 lnp 2 , and 
for functions F going like (\n\np 2 ) N /p 2 , K>F behaves like (\n\np 2 ) N+1 /Np 2 \np 2 . One can then define the inverse 
(1 — Ky)^ 1 by the pedestrian means of summing the series in Eq. (|4ip . with the result (valid in the UV): 

Q(P 2 ) = explaining] - ^^[lnp 2 ] Q . (45) 

pZ | n pZ pi m pA 



It is clear that the meaning of lnlnp in the UV is 



in view of: 



lnlnp 2 ^ln[^-] (46) 
9 (P ) 



= 1 + bg 2 ln[l + -^] (47) 



g 2 {p 2 ) Am 2 

which follows from Eqs. (I17|) . (|17p . Then the result for the UV behavior coming from combining confinement and 
one-gluon terms is: 

This UV behavior is what the RG dictates, and by adding in confinement effects we are able to give the prefactor al c . 
From the work of [23[ this allows an estimate of the (qq) condensate, although we will not pursue that further here 
because of various complications (see, for example, [26|]). Finally, we can change M c (p 2 ) to Mi c (p 2 ), which is M c (p 2 ) 
modified by the IR one-gluon corrections (c/. Eq. (134)) to incorporate the IR corrections from one-gluon exchange. 
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VIII. THE INTEGRABILITY CONDITION 



Is it possible to impose an integrability condition on a confining effective propagator, of the form: 

J d 4 fcD e //(p -k) = const. (49) 

where the constant is finite and independent of any regulator masses that were introduced? If so, the explicit regulator 
m can be dropped, and the integrability condition provides an implicit regulator that makes the Minkowski-space 
JBW equation finite. 

The criterion of integrability certainly fails for the Feynman propagator l/[(p — fc) 2 — m 2 + it] 2 , whether or not 
there is a regulator m. But it does hold for the principal-part propagator. The principal part propagator is physically 
interesting for other reasons, because it is real and leads to no decay processes for the supposed area of the area 
law. It is the sum of the advanced and retarded propagators, discussed long ago by Wheeler and Feynman in their 
preliminary formulation of QED. 

Begin with the standard principal-part propagator known as D(x) with mass m, defined as: 

D(x) = l[D adv (x) + D ret (x)} = -j^-rj ( d 4 ke- lk - x P 1? 1 • (50) 



D{x) = —5{x 2 ) + ^^pij^mVx^). (51) 

47T j? ^ 



in space-time it has the value: 



To go to a confining potential, multiply by 8nKp and differentiate with respect to m 2 



D eff (x) = -8*^ (52) 



<TkP 



-8ttK f 
(2tt) 4 J " [{k 

~ KF ' "<r 2 \ T„(m^f^2\ 



2 



g -ifc, 



(x ) J (mVx 2 ). 

The static potential, the time integral of D e ff, is the usual Kpr plus an infinite constant coming from the integral 
over all time. So we have not entirely exorcised the singularities of a confining potential. Earlier it was argued 
that this singular term can be gauged away, with a singular gauge transformation. 

The integrability condition will be studied in detail in a later work. For now it is enough to say that if in the 
Minkowski-space gap equation we make the drastic approximation 9(x 2 ) — > 1 and certain other simplifications, the 
result is a wrong-sign Ai 4 theory with a homogeneous solution M. 2 — Kp/ir. Here M(x) is the Fourier transform of 
M(p 2 ). 

A potential problem shows up: How to continue such a propagator to Euclidean spacetime for use in a Euclidean 
JBW equation? There is no obvious continuation of the theta-function. If there is a reasonable continuation, will it 
preserve an integrability condition? We have found infinite classes of (the S-wave projections of) Euclidean propagators 
with no regulator mass that satisfy an integrability condition and behave as 8nKp / k for certain ranges of momenta 

1/2 

< Kp , and found that the integrability condition insures singularity-free CSB. One example: In the integral equation 



(l27t multiply the integrand of J<(p ) by p /k , and then set m = in both J< and J>. The resulting massless kernel 
yields a finite constant when integrated over k 2 dk 2 , and the corresponding gap equation has properties very similar 
to those found in connection with the original kernel of Eq. (I2T1) . 



IX. SUMMARY 



We suggest that when chiral symmetry is broken, leading to a running quark mass M(p 2 ), the ensuing massless 
Goldstone bosons contribute to an entropy-driven condensate such as (qq) by ramifying a large number of branches 
from a basic "trunk" Wilson loop that itself shows large fluctuations. Even if the Wilson loop represents quenched 
quarks, and so is incapable of breaking, configurations with the q and the q far apart are quite improbable compared 

— 1/2 

to ramified configurations where they are only separated by a distance of order K F , because the associated area- 
law action is large compared to the entropy. This means that in a linearly-rising potential Kpr, separations with 
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r 3> K F are not probed, and the potential at such large distances is irrelevant. In turn, this means that confinement 

— 1/2 

dynamics need a linearly-rising potential only out to a finite distance r ~ K F 

From these considerations we formulate a non-singular confining gap equation of JBW type, with a mass m in the 
confining effective propagator treated not as a regulator to be set to zero, but as a finite mass ~ M(0) that can be 
estimated from processes that respect the Abelian gauge invariance associated with the confining effective propagator. 
The static potential coming from the confining effective propagator rises linearly only out to a finite distance, and 
has a negative term at the origin that we identify with entropic effects. We estimate m through comparison with an 
extension of an old calculation (having Abelian gauge invariance) in which m could be properly treated as a regulator 
mass, to be sent to zero after cancellations. The extended calculation replaces the regulator by a specific and physical 
mass, whose dynamical effects are equivalent to keeping m as finite and of this value in the JBW equation. 

We also studied one-gluon effects within the same framework, thereby finding the correct large-momentum behavior 
of QCD as known from the RG, but with a calculable prefactor. The final result, including IR one-gluon enhancements, 
is a quark mass M(0) « 250 MeV. This is somewhat smaller than the commonly-quoted value of 300 MeV, which is 
largely based on not on true dynamical estimates but on assuming that the sum of quark masses represents most of 
the mass of the hadron in question. Our mass has a different interpretation, given its entropic underpinnings. 

Much remains to be done, both in formulating and solving more elaborate (and more accurate) forms of the gap 
equation, including in Minkowski space; relating these to pion dynamics; and perhaps making some progress in 
understanding entropic effects quantitatively. 
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